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Abstract 


Microreactor technology is a promising approach in harnessing the high energy density of hydrocarbons and is being used to produce hydrogen- 
rich gases by reforming of methanol and other liquid hydrocarbons. However, on-demand H, generation for miniature proton exchange membrane 
fuel cell (PEMFC) systems has been a bottleneck problem, which has limited the development and demonstration of the PEMFC for high- 
performance portable power. A number of crucial challenges exist for the realization of practical portable fuel processors. Among these, the 
management of heat in a compact format is perhaps the most crucial challenge for portable fuel processors. In this study, a silicon microreactor- 
based catalytic methanol steam reforming reactor was designed, fabricated, and demonstrated in the context of complete thermal integration to 
understand this critical issue and develop a knowledge base required to rationally design and integrate the microchemical components of a fuel 
processor. Detailed thermal and reaction experiments were carried out to demonstrate the potential of microreactor-based on-demand Hp generation. 
Based on thermal characterization experiments, the heat loss mechanisms and effective convective heat coefficients from the planar microreactor 
structure were determined and suggestions were made for scale up and implementation of packaging schemes to reduce different modes of heat 


losses. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


As fuel cell (FC) gain attention as potential energy solution 
of the future, one of the most promising fields of application 
is portable power, as a prominent alternative to batteries [1]. 
Batteries are a multi-billion dollar market and the technology 
has made significant strides in the last couple of decades to 
meet growing power demand of advanced portable computing, 
handheld personal devices, and other microelectronics systems. 
The rechargeable battery technology has progressed from nickel 
cadmium (NiCD) to nickel meta hydride (NiMH) to lithium 
ion (Li-Ion). Until recently, the gradual improvement of Li- 
based batteries was enough to power mobile computing devices. 
Unfortunately, the energy density available from current bat- 
tery technology is inadequate to match the power demands 
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of emerging portable devices [2]. This need has encouraged 
researchers to continue to vigorously explore FCs as miniature 
power sources. 

High energy density is required for advanced portable appli- 
cations for consumer use (e.g., laptop computers, PDAs, cellular 
phones, medical and telecommunication devices), and for mili- 
tary and intelligence use (e.g., remote sensors (<20 W), systems 
for the soldier (<500 W), silent power generation and battery 
recharging (200 W-2kW), and mobile power generation and 
auxiliary power units (>2kW)) [3], with the goal of allowing 
these devices to operate for longer times with less recharging. 
Miniatures FCs offer promise for the conversion of chemical 
energy into electrical energy in a highly efficient way achiev- 
ing the goals of power dense, light weight, and compact power 
sources. 

The proton exchange membrane fuel cell (PEMFC) is partic- 
ularly attractive and promising for portable applications because 
of its ability to provide high power density and nearly instant 
power. Moreover, its simplicity in design and operation allow- 
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Nomenclature 

A surface area (m?) 

Aw cross-sectional area of the wire (m?) 

CMeoH methanol concentration in feed (mol cm~?) 

dy hydraulic diameter (m) 

Dap diffusivity of methanol in gas mixture (m? s~!) 

E Emissive power (W m~?) 

Ep. spectral blackbody emissive power 
(Wm? pm!) 

h heat transfer coefficient (W m~? K7!) 

I current (A) 

Å, radiative intensity (W m~? s7!) 

kg mass transfer coefficient (m s7!) 

KY spectral absorption coefficient of Pyrex™ (m7!) 

ke spectral absorption coefficient of silicon (m7!) 

l length of the wire (m) 

Nyu Nusselt number 

NRa Rayleigh number 

P' perimeter of the wire (m) 

Qconv. heat loss by convection 

Qin heat input (W) 

Qloss total heat loss (W) 

Qrad. heat loss by radiation 

Qwires heat loss from the attached wires 

TMeOH methanol reaction rate (mol cm~? s7 1) 

Sh Sherwood number 

Ts surface temperature (K) 

Tx ambient temperature (K) 

V voltage (V) 

Greek letter 

a constant used in equation (5) 


ing fabrication of a compact and lightweight cell, mild operating 
conditions (temperature less then 100°C, pressure in the range 
of a few atmospheres), and relatively low loading of precious 
metal catalysts are additional strong advantages [1,4]. 

The success of PEMFC technology for portable power 
depends heavily on the development of efficient ways of deliv- 
ering the appropriate fuel to the cell. The PEMFC operating 
on hydrogen as the fuel holds the most promise, offering the 
best possible efficiency and power density [5,6]. However, while 
hydrogen offers high energy density on a per weight basis, it has 
a poor energy density per volume. Gaseous hydrogen is there- 
fore difficult and inefficient to store and several high energy 
density storage methods have been investigated such as highly 
pressurized cylinders, metal hydrides, chemical hydrides, and 
even carbon nanotubes [7,8]. Since on-board storage of hydro- 
gen is unsafe and energy inefficient, two main alternatives are 
commonly envisioned and employed for portability. 

One approach, the direct FC employs fuels like methanol, 
ethanol, formic acid, etc., directly in the cell to generate elec- 
tricity. The direct FC has the advantage of a smaller total system 
and weight, as it does not need an in-situ fuel processing (FP) 


unit to extract hydrogen. Among all direct FCs investigated 
at the portable scale, more research and development efforts 
are focused on the development of direct methanol fuel cell 
(DMFC). Several groups are pursuing low power DMFCs for 
portable power applications [9]. Although the DMFC offers 
promising possibilities as a miniature power source, the technol- 
ogy has to overcome a number of challenges to be competitive 
with battery technology in terms of efficiency, power density, 
and cost. Methanol crossover and the relatively slow reaction 
rate of methanol oxidation to hydrogen ions and carbon diox- 
ide at the anode are persisting challenges that currently limit 
the DMFC to less stringent applications. In addition to DMFC, 
the direct ethanol fuel cell (DEFC) and direct formic acid fuel 
cell (DFAFC) are also being developed. Ethanol is considered 
an attractive fuel mainly because of its higher theoretical energy 
density than methanol, non-toxicity, availability, and renewabil- 
ity [10,11]. However, the available power density from DEFC 
is currently about half that of DMFC resulting from relatively 
slower electrooxidation reaction kinetics at the anode. Recently, 
the DFAFC has been demonstrated as an attractive alternative 
[12]. The theoretical energy density of formic acid is less than 
one-third of methanol. However, formic acid has other advan- 
tages like faster anode kinetics and smaller crossover through 
Nafion membranes. Unlike DMFC where methanol crossover 
limits the use of higher methanol concentration solutions, it is 
possible to run the DFAFC with high concentrations of formic 
acid. However, it was recently observed that the performance of 
the cell dropped as the formic acid concentration was increased 
beyond 10M because of decreases in anode catalytic activity 
and an increase of the fuel crossover [13]. 

The other, higher performance alternative is to generate pure 
hydrogen on site and on demand by processing easily stored, 
high energy density liquid fuels [14]. FP adds complexity, 
weight, volume, and cost compared to the DMFC or a sys- 
tem running directly on hydrogen, however the effective energy 
storage density is significantly higher. 

A wide variety of energy dense hydrocarbons such as nat- 
ural gas, methanol, gasoline, and diesel can be reformed into 
hydrogen-rich streams for FCs. Different processes are available 
for hydrogen extraction and subsequent purification, there- 
fore, the type of hydrocarbon and method used to produce the 
hydrogen play a major role in the design of a FP system. Pro- 
cessing of liquid hydrocarbons to hydrogen involves a series 
of steps including fuel vaporization, reforming (three main pro- 
cesses available are partial oxidation, steam reforming (SR), and 
autothermal reforming), water gas shift (WGS), and CO clean up 
(preferential oxidation (PrOx), hydrogen separation by perme- 
ation through a hydrogen-selective membrane, or pressure swing 
adsorption) before feeding hydrogen rich gases to the FC. 

Methanol is the fuel of choice for mobile portable power 
application as it is sulfur free and can be reformed at rela- 
tively low temperatures (250-300 °C), simplifying the micro 
FP design. Moreover, it offers other advantages like high 
hydrogen—carbon ratio, high energy density, and ready avail- 
ability. 

The following chemical reactions take place during the SR 
process. 
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CH30H + H20 + CO2+3H2 (AHR= 50kJ/mol) (1) 


CH30H < CO + 2H2 (AHR= 90kJ/mol) (2) 


CO + H20 © CO2+H2 (AHR= —41kJ/mol) (3) 


Eq. (1) is the main reforming reaction, Eq. (2) is the methanol 
decomposition reaction, while the third equation represents the 
water—gas shift reaction. The reforming reaction is endother- 
mic in nature and thus requires external energy. Also energy 
is required to vaporize and superheat the methanol—water mix- 
ture to the desired SR reaction temperature. These requirements 
can be reasonably met by combusting excess hydrogen from the 
anode exhaust of the FC. The catalytic SR of methanol on a 
CuO/ZnO based catalyst is a highly efficient process, producing 
a reformate with high hydrogen concentration (65-70%) while 
maintaining high CO) selectivity compared to other available 
reforming processes. A CO concentration of 1-2% is attained 
while operating at temperatures ranging from 220 to 260°C as 
a consequence of reaction thermodynamics [15,16]. Therefore 
a methanol-based FP system must include a CO clean up step, 
as this relatively low level of CO is still enough to poison the 
PEMFC catalyst. However, the low CO concentration eliminates 
the need for WGS, requiring only a PrOx step to reduce CO 
below 10 ppm. 


2. Microscale FP: issues and challenges 


A block diagram for an overall system was developed and 
is shown in Fig. 1. The system burns the unconverted hydro- 
gen from the FC anode in the combustor to supply the heat 
required for vaporizing the methanol-water feed and to carry out 
the endothermic reforming reaction. The total system includes 
key microreactor components (combustor, vaporizer, SR, and 
PrOx reactor), FC stack, a diluted methanol storage tank (water 
recycling is not employed for simplicity), a rechargeable battery 


H,, CO,, CO, H,O 


Steam reformer 


for start up and rapid load following, insulation, and peripheral 
components. In addition, a separate storage of pure methanol 
is required as an additional fuel in the combustor if the uncon- 
verted anode hydrogen does not meet the energy requirement 
for vaporizer and SR. 

Considerable efforts have been made to develop an integrated 
FP/FC system, however, several issues and challenges persist. 
Based on the literature and our own research, many of these 
key issues are identified and discussed in our previous paper 
[17]. Among them, miniaturization of system components and 
thermal management are perhaps the most crucial. In order to 
understand the relationships and tradeoffs of these issues, a study 
encompassing energy density, component sizing, and thermal 
integration of key FP components was performed. 


2.1. Analysis of energy density with methanol-based FP/FC 
system 


Energy density is the ratio of energy available from a power 
source to its volume (usually quoted in Wh17!) or weight 
(Whkg~'), while for a fuel it is the amount of energy that can 
be contained in a specific quantity of the fuel itself. 

Pure methanol has an energy density of 4380Wh17! 
(5600 Whkg~!), which is about 15 times the projected energy 
density of a Li-ion secondary battery (300 Wh17!), meaning 
that the FP/FC unit has superior energy density even with an 
overall conversion efficiency as low as 7%. In the case where 
no water recycling is employed to minimize system complexity, 
water has to be carried with methanol for the reforming. With a 
stoichiometric steam to carbon ratio of 1:1, this reduces the net 
energy density to 3000 Wh17! (4550 Wh kg7!). 

The energy density comparison between methanol-based FC 
technology with batteries illustrates the potential of this tech- 
nology for the consumer portable market. If an advanced cell 
phone device consuming 1 W is powered by a 50 ml methanol 


Recycled H, 
= ro | Lt 


Fuel cell stack 


Electrical 
energy 


Combustor 


Methanol, water 


storage tank 


Fig. 1. Block diagram for converting methanol to electricity using an integrated FP/FC system. 
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Table 1 

Energy density at various mission durations 

Mission duration DARPA goal Weight (kg) of FP/FC* Weight of BA Weight ratio 
(Whkg~!) (energy density, Wh kg~!) 5590? (kg) (FP/FC)/BA 5590 

3h, 60 Wh 1000 0.55 (110) 0.35 1.57 

3 days, 1440 Wh 2000 1.75 (820) 8.32 0.21 

10 days, 4800 Wh 3000 4.70 (1020) 27.75 0.17 


Based on DARPA expectations and estimated from methanol-based FP/FC integrated system. 


2 Weight of FP/FC excluding (methanol + water) storage: 500 g. 


P BA 5590, a lithium-ion battery (energy density of 170 Whkg~!) is a widely used battery for military applications. 


cartridge, and considering only 10% overall efficiency, a con- 
tinuous transmission time of 22 h without recharge is possible. 
This represents a significant gain relative to current capabili- 
ties. 

A number of military applications have been identified 
which represent a significant market for miniature power. The 
Department of Defense (DoD) seeks lighter and more compact 
electrical power sources for individual soldier, robotic, and other 
combat missions. The higher energy density will allow soldiers 
to carry less weight and to extend mission duration. 

In 2002, DARPA set the energy density goals for three mis- 
sion scenarios at an average power level of 20 W [18], which are 
listed in Table 1. The expected specific energy includes the total 
system and fuel. 

If we assume an optimistic conversion efficiency (which 
includes FP conversion efficiency and FC efficiency) of 25%, 
available energy density would be 1140 Whkg™!. It is pos- 
sible to determine whether methanol-based systems can meet 
energy density expectations set by DARPA for different mis- 
sions assuming the conditions listed in the caption to Table 1. 
This comparison is shown in Table |. We see clearly that DARPA 
has set aggressive energy density expectations for these mis- 
sions, which can only be achieved with diesel or other high 
energy density hydrocarbons; or perhaps with alternative tech- 
nologies. 

Table 1 also shows the comparison between the weights of 
the BA 5590 battery versus a methanol FP/PEMFC system. 
Net energy density increases favorably for extended missions 
since only fuel is added to increase the duration. As shown in 
Table 1 missions of 10-day duration result in a net energy den- 
sity of 1020 Whkg™!, which is a six times improvement over 
the BA5590 (170 Whkg7!). 


2.2. Miniaturization: microchemical systems for microscale 
FP 


Microchemical systems have the potential for highly com- 
pact power sources through integration of a FC with all the unit 
operation components of a FP along with microstructured sen- 
sors, actuators and other “balance of plant” (BOP) devices. [19]. 
Their microscale dimensions result in high heat and mass trans- 
fer rates [20], making them efficient in applications where it 
is necessary to thermally couple endothermic and exothermic 
processes. The excellent transport and low thermal mass ensure 
short times to reach working temperature at start-up. 


2.3. Miniaturization: component sizing and compatible 
catalyst technology 


The initial objective of our group’s research efforts was 
to perform a model-based design and integration study of 
microchannel reactor components of a fuel-processing unit for 
converting methanol to hydrogen with <10 ppm CO for a 20W; 
PEMFC. The ASPEN PLUS™ process simulator was used to 
simulate the overall FP/PEMFC system. The ASPEN analysis 
provided major operating parameters such as temperatures, pres- 
sures, and flow rates needed to carry out a preliminary design 
for key FP microreactor components as a first order estimate of 
the size requirements for the 20W. micro FP [21]. 

The SR, combustor, and PrOx microreactors were designed 
based on a thin-film catalyst approach [22], which consists of 
a ~10 pum layer of catalyst deposited on microchannel walls. 
For the preliminary design of the SR and combustor, a kinetic 
model based on power-law rate expressions from the literature 
were used [23,24]. The PrOx microreactor design is based on 
our own kinetic experimental data as our group has successfully 
demonstrated a 0.25W, single channel silicon microreactor with 
Pt/Al 03 thin-film wall catalyst for PrOx of CO in hydrogen rich 
reformate [25]. This modeling approach estimated the catalyst 
loading for each of the microreactor components, which helped 
determine other microreactor design details such as size and 
number of microchannels required, overall reactor layer size 
and geometry, and the number of layers to be stacked, etc. The 
process parameters and the design results for each component 
are listed in Table 2. 

Miniaturization of each component by reducing weight and 
volume and high-energy efficiency were the main objective of 
this work. Due to different demands placed on the different 
reaction steps, there are different requirements for the catalysts 
for each of the components. The sizing estimates show that the 
thin film catalyst approach can be used for PrOx and combustor 
microreactors but it is not favorable for the SR microreactor at 
this scale, based on the achievable loading. Note that the com- 
bustor required the least amount of catalyst loading, such that a 
single Pt wire would be sufficient. 

The results obtained from the above 1D sizing-analysis 
showed that the SR is the largest component and it drives the 
overall size of the micro FP. The SR reaction is relatively slow 
and hence needs larger amounts of catalyst at this temperature 
range to achieve high conversion. It is not possible to load a suf- 
ficient catalyst amount to further miniaturize the SR if the thin 
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Table 2 
Preliminary sizing estimates: process parameters and design results 


Component Process parameters and design results 


SR We, 60% FC efficiency, 75% H3 utilization 
MeOH conversion 100% 
Fyeou = 6.233e—5 mol s7! 
Fu,0 = 8.10e—5 mol s7! 
H20:MeOH = 1.3:1 
260°C 
latm 
Microchannel: 50 wm x 400 pm 
Parallel microchannels: 1000 
Channel to channel spacing: 50 ym 
Reaction zone width: 10 cm 
Catalyst film: Cu/ZnO/Al203, 5 pm thickness 
(0.0478 g cm7!) 
Reactor layers: 24 
Total catalyst loading: 2.87 g 


PrOx 160°C 
Catalyst film: Pt/Al203, 5 pm thickness 
(0.0478 g cm7!) 
CO in reformate: 2.2% 
Reactor layers: 9 
Total reactor loading: 1.1 g 


Off gas recycle to combustor: 2.5e—4 mol s7! 
Gas composition: H2: 33.8%, CO2: 38.8%, 
H20: 11.8%, O2: 0.6%, and N2: 15.1% 

Total catalyst loading: 1.2 mg 


Combustor 


film approach is used. Thus a different approach to incorporate 
catalyst into the small reaction zone of microreactor at higher 
loading is needed. Packing the reformer with catalyst powder 
and use of a foam-based cellular catalyst [26], offering signifi- 
cantly higher loading than thin-films are other options that were 
also considered. 


2.4. Thermal management in microscale FP 


Thermal management in miniature systems is perhaps the 
most crucial challenge for microscale FPs. Thermal coupling 
is needed to allow transfer of energy from the heat produc- 
ing combustor to the endothermic SR. Coupling endothermic 
and exothermic components and minimizing losses promotes a 
high thermal efficiency. However, such coupling must be accom- 
plished in a manner that permits the maintenance of specific 
temperatures in the various components as shown in Fig. 2 and 
maintains the surface of the package near room temperature. 
Microreactors generally offer high heat transfer rates mainly 
because of high surface-to-volume ratio, short conduction paths, 
and high thermal conductivity of materials used. This charac- 
teristic results in efficient heat transfer between the different 
processing steps but at the same time results in higher heat 
losses to the ambient. Therefore, thermal management offers 
a dual challenge of opposing the heat losses from the system 
that arise from high surface-to-volume ratio in conjunction with 
maintaining temperature gradients within the system to allow 
desired conditions in the unit reaction steps. 

Different materials have been used for the construction of 
microscale FP components, including metal [27-29], silicon 


H, CO, 
PrOx, 155°C 


Methanol/ water 
SR, 260°C 


Combustor, 350°C 


Air Recycled H, 


Fig. 2. Integration of fuel processing components. 


[30-33], ceramic [34], and glass [35]. However, all the materials 
approaches are confronted with the thermal management issues 
described above. Because microfabrication processes leverage 
two-dimensional replication technologies like photolithography, 
a flat, sheet-like architecture is the norm for microfabricated 
components of FPs. This kind of low aspect ratio, planar archi- 
tecture offers an excellent inter-component heat transfer, but is 
poor for heat retention and results in significant losses to the 
ambient. 


3. Thermal management: demonstration of an 
integrated SR 


To address the crucial issue of heat management, we present 
the development of a silicon microreactor based methanol SR 
with emphasize on thermal integration of the overall system. As 
shown in Fig. 3, a silicon micro SR was designed and integrated 
with microfabricated insulation as dictated by an overall system 
design. An integrated design was made where thin metal film 
heaters and temperature sensors were directly embedded into 
the unit. External heat was introduced to simulate integration 
with a combustor. Detailed thermal data (heat input required, 
heat losses, temperature profiles) and reaction data (conversion 
and selectivity) were acquired. 

Because of its various advantages, silicon microfabrication 
technology was employed in this work to achieve minia- 
turization. However, our intent was to perform thermal and 
reaction experiments whose outcome depends most critically on 
microscale geometry and less so on the materials of construction. 
In this way the results obtained may be generalizable to any of the 
materials technologies. Silicon is a reasonable material choice 
as it is lighter than stainless steel and possesses high mechanical 


155°C 


Microfabricated Insulation 
Microfabricated Insulation 


Thin-film heater (350°C) 


Methanol/water T}, CO, CO, 


Fig. 3. Demonstration of thermally integrated SR. 
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strength, chemical inertness, thermal stability, and large thermal 
conductivity. With the established silicon microfabrication tech- 
nology, it is possible to create novel microchannel and reactor 
configurations in the sub-millimeter range. An important advan- 
tage of silicon microfabrication is its compatibility with thin 
film technology. This facilitates eventual on-chip integration of 
functional elements like sensors and actuators allowing compact 
process control subsystems. 


3.1. Design of a micro methanol SR 


A silicon microreactor with a reaction zone of 1 cm x 1 cm x 
400 um (0.04cm?) was designed. The design incorporates a 
packed bed of commercial copper-based powder catalyst. The 
packing of the reaction bed resulted in a catalyst loading of 50 mg 
calculated using the bulk density of the catalyst. A kinetic model 
was formulated using an Arrhenius type power-law rate expres- 
sion for methanol reforming over a Cu/ZnO/Al203 catalyst 
[23], assuming no mass transfer resistance. Using this model, 
with the knowledge of reaction temperature and catalyst load- 


(a) Flow manifold 


Methanol/water 


ing, methanol and water flows (S/C molar ratio of 1.3:1) were 
determined to achieve essentially 100% conversion of methanol. 
Using a methanol flow rate of 1.32e—4 mol min! at 260 °C, 
the reaction produces 3.84e—4 mol min™! (9 sccm) of hydro- 
gen, which was sufficient to produce 0.65 W from a FC assuming 
60% efficiency and 75% hydrogen utilization in the FC. 

Fig. 4(a) shows the design of the SR microreactor. The reac- 
tor consisted of an inlet to introduce the mixture of vaporized 
methanol and steam, an outlet to remove reformate, a packed 
bed reaction zone, flow manifold structure, and filter structures 
at the outlet to trap catalyst particles. 

Pressure drop for the reaction zone packed with catalyst was 
estimated from the Ergun equation [36]. For catalyst particles 
of 75 pm, the reactor of 1-cm bed length resulted in a nearly 
insignificant pressure drop of 190Pa. This low pressure drop 
is consistent with the low superficial velocity resulting from the 
flow rates of gaseous reactants through the 1-cm wide x 400 um 
high cross-section. Low pressure drop across the bed can lead to 
maldistribution of flow and the resulting non-uniform catalyst 
utilization and residence time. In order to avoid this situation, 


Catalyst traps 


Hp, CO, CO» 


KS/RSB 


Reaction zone (catalyst packed bed) 


(b) 


Fig. 4. (a) Design of the SR microreactor. (b) 2D axisymmetric flow model used in CFD simulation and velocity distribution in microreactor. 
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an inlet manifold structure was designed to obtain uniform flow 
distribution across the reaction zone width. The manifold chan- 
nels were designed to ensure that catalyst particles easily flow 
into the reaction zone during loading. The FLOTRAN CFD sim- 
ulation package in ANSYS Multiphysics software was used to 
determine the velocity profile across the reaction zone for a given 
boundary condition of flow velocity at the inlet with an outlet 
boundary condition of atmospheric pressure. Fig. 4(b) shows the 
2D axisymmetric flow model and the velocity profile along the 
packed bed. 


3.2. Internal heat management: thermal coupling and 
design of insulation 


In order to maintain desired temperature gradients between 
FP components, insulators are required to bridge the temper- 
ature differences. The goal was to design insulation such that 
it adds little to the weight and volume of the total system and 
can be integrated with microfabrication. The high surface-to- 
volume ratio of microreactors, high thermal conductivity of 
silicon, and sub-millimeter conduction paths make it difficult 
to obtain the desired thermal isolation between components. It 
is therefore, necessary to employ special structures like thin sus- 
pended membranes, vacuum gaps, or insulating materials with 
ultra low thermal conductivity to bridge these temperature dif- 
ferences and to minimize conductive and convective heat losses 
to the surroundings. 


(a) 1.000 


0.800 


The best commercial insulators have thermal conductivities 
in the range of 0.02 W mK~! and are difficult to integrate into a 
microfabricated unit. However, vacuum packaging of microre- 
actors can provide an effective means of insulation. From the 
kinetic theory of gases, it is known that the thermal conduc- 
tivity of a gas is approximately independent of pressure for 
atmosphere and above [37], and decreases with sub-atmospheric 
pressures, as the mean free path becomes less than the enclo- 
sure dimensions [38]. Fig. 5(a) shows a plot of kp/ko (kp being 
thermal conductivity of air at sub-atmospheric pressure [39] 
and ko being thermal conductivity of air at atmospheric pres- 
sure) versus sub-atmospheric pressures for a cavity depth of 
500 um. A microfabricated cavity surrounding the reactor filled 
with low-pressure gas offers low thermal conductivity down to 
approximately 0.001 W mK~!. Silicon microfabrication enables 
a straightforward approach to this structure, as the cavity etched 
in silicon can be sealed to trap the desired level of vacuum and 
the desired residual gas species. Fig. 5(b) shows the concept for 
integrating vacuum insulation. 

3D thermal simulation using ANSYS® simulation software 
was carried out to design microfabricated vacuum insulation 
such that it can maintain the required temperature gradients. 
The simulation also determined placement of the combustor and 
PrOx zones in reference to the reformer reaction zone, pres- 
sure in the cavity and the depth of cavity required. As noted 
earlier, the reforming microreactor was designed to produce 
3.84e—4 mol min! of hydrogen. Considering 75% hydrogen 


0.600 
l d = 0.5 mm 
ky (air at 760 torr and 533 K) = 0.04 W/mK 
0.400 
0.200 
0.000 
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(b) Vacuum cavity, P = 0.001-0.1 torr Glass (Pyrex™) 
SR reaction 
zone 


0.5 mm mm 


if i 


a 


Fig. 5. (a) Thermal conductivity of air at sub-atmospheric pressures for a cavity depth of 500 um. (b) Vacuum insulation chip with a single microfabricated cavity 


in Si, and reformer reaction zone sandwiched between insulation chips. 
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(a) 155°C 


Methanol/steam 
— > 
(m.Cp-Tin) = 0.038W 


Microfabricated Insulation 


SR (-0.13W) 
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(b) Combustor Zone 


Microfabricated Insulation 


Í 0.38W 
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210.994 


Fig. 6. (a) The concept of 3D thermal model used to design microfabricated vacuum insulation. (b) 3D thermal simulation results showing the temperature distribution 


in different compartments. 


utilization in the FC, the remaining 9.6e—5 mol min”! of hydro- 
gen in the anode effluent when burned, gives 0.38 W based on 
the lower heating value of hydrogen. The temperature bound- 
ary conditions, endothermic heat of reforming, and reactant 
and product sensible heats used in the 3D model are shown 
in Fig. 6(a) while Fig. 6(b) shows the result. The simulation 
yielded the temperature distributions in each of the compart- 
ments. Thermal conductivity of air in the cavity was varied 
until the required component temperatures were achieved. For 
generality, the model did not consider heat losses. 


3.3. Design and fabrication of an integrated SR 


The design of the integrated SR is shown in Fig. 7(a). The 
stack of six chips consisted of insulation chips surrounding 
the reforming microreactor. The silicon reactor was covered on 
the top with a Pyrex™ sheet of the same lateral dimensions. 
Thin film resistance temperature detector (RTD) sensors were 
patterned at different positions on the backside of the silicon 
microreactor. Three sensors were placed at different positions 
to measure the reaction zone temperature distribution and one 
each was placed at the inlet and outlet section as shown in the 
superimposed image (Fig. 7(b)) of the front and backside CAD 
masks of the microreactor. The sensors have meandering resis- 
tors with a four-wire design to minimize lead resistance effects. 
The insulation consisted of an assembly of silicon and glass 
chips. The silicon chip had a cavity in the center and was bonded 
with Pyrex™ heater and temperature sensor substrate as shown 


in Fig. 7(c). The relative sizes of the combustor and PrOx units 
were determined from the initial size estimation discussed in 
Section 2.3. The insulation requirements, described next, were 
determined based on thermal and structural considerations. 

The insulation requirement between the combustor and the 
reformer is not demanding, as the reformer needs endother- 
mic heat of reaction. An appropriate thermal conductivity and 
insulating thickness is obtained in this case by a 50-um deep 
cavity containing air at atmospheric pressure. The lateral dimen- 
sions of the silicon and glass are thus 1.1 cm x 1.4m while the 
cavity is slightly smaller at 1.0cm x 1.3 cm. The combustor is 
2.8 mm x 2.8 mm and is centered over an insulation chip keeping 
a lateral distance of 7.2mm x 10.2 mm from conducting verti- 
cal paths. An integrated on-chip thin-film Pt heater meandering 
along this area is used to simulate the combustor. Three thin- 
film temperature sensors are also patterned on the glass chip. 
One sensor was located right in the center of the combustor 
zone, while two others were placed on either side at a distance 
of 600 um from the center. 

On the opposite side of the assembly, more stringent thermal 
isolation is required between the reformer and the PrOx. For this 
reason, the insulation between the reformer and the PrOx had a 
200-um deep cavity and was sealed with vacuum of less than 
5.107 Torr. The silicon insulation chip of lateral dimensions 
1.2cm x 1.5cm had a cavity of 1.1 cm x 1.4cm. The effective 
insulation in this case requires higher vacuum, a deeper cav- 
ity, and smaller conductive distance between the cavity and the 
insulation chip. Since silicon has a high thermal conductivity 
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(a) Insulation 
200-um depth, P = 0.005 torr 


k = 0.0006 W/mK T Sensors 500 
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Pyrex™ Outlet 
Inlet -= 
Heat Insulation 
CALEL 50-um depth, Atm. pressure 
k = 0.04 W/mK 


(b) 


(c) 


Fig. 7. (a) The design of the integrated SR. (b) Superimposed image of the front and backside CAD masks of the microreactor. (c) Mask design for silicon and glass 
chips of the insulation layer. Silicon chip with cavity to be etched in the center and glass chip patterned with thin film heater and temperature sensors. 


(150 W mK~!), only a 500-um wide border was maintained 
around the cavity on the silicon insulation chip to minimize 
vertical conduction of heat. 

The effects of pressure imbalance across relatively thin slabs 
must be considered as this situation can result in stress accumu- 
lation and deformation of silicon or Pyrex™. In extreme cases, 
fracture can result during wafer dicing to separate individual 
devices or later during device operation. Structural simulation 
was carried out with ANSYS® to determine the deformation 
and the stress formation in silicon and glass for the required 
vacuum inside the cavity. This simulation indicated that an 800- 
um thick wafer with 200-um cavity depth would suffice for 
PrOx insulation. For the designed insulation cavity chip size 
of 1.1 cm x 1.4cm, the pressure force resulted in a predicted 
deformation of less than a micron in silicon and glass. 

Four microreactors were fabricated on a single silicon wafer 
(p-type <100> 4-in. diameter, 550-~m thick, double-side 
polished) by silicon bulk micromachining techniques. Pho- 
tolithography was employed first to pattern the front side of 
the microreactors using thick photoresist. These structures were 


then etched to a depth of 400 um using deep reactive ion etch- 
ing (DRIE) by inductively coupled plasma (ICP) [40]. A 2-um 
thick silicon dioxide hard mask was used for deep silicon etch- 
ing. Once the front side was processed, a second mask was used 
for backside patterning of inlet and outlet vias. These vias were 
then opened from the backside by wafer through etching with 
DRIE. Thin-film Pt temperature sensors were then patterned on 
the backside of wafer. A 10 nm Ti adhesion layer was deposited 
followed by 200 nm of Pt without breaking the vacuum of e- 
beam evaporation. The photoresist was dissolved, removing the 
excess Ti/Pt by lift-off, leaving the temperature sensors. The 
fabricated silicon wafer was then sealed with Pyrex™ glass by 
anodic bonding. The bonded reactor wafers were then diced into 
individual microreactors. 

The vacuum insulation was formed by creating first a square 
cavity pattern on the silicon wafer and etching to a required depth 
using DRIE. Thin-film heater and temperature sensors were fab- 
ricated on the Pyrex™ glass wafers using the lift-off technique. 
Here, the metal lines consist of a layer of Ti (10 nm) and Pt 
(300 nm). The patterned glass wafer was then annealed at 650°C 
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zone. 


Temp. sensor (Te) in the 
center of the combustor 


Packed bed of catalyst 


SEM of thin film Ti/Pt 
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Fig. 8. (a) Frontside of silicon SR microreactor with SEM images of inlet via, flow manifold, and catalyst traps. (b) Backside of SR microreactor showing patterned 
temperature sensors. (c) Incorporation of catalyst in the form of packed bed by vacuum loading. Catalyst loading achieved was 53.1 mg. (d) PrOx side insulation 
chips. Etched cavity in silicon bonded under vacuum with patterned glass. (e) Combustor side insulation chips with SEM images of central combustor (heater) zone 


and temperature sensors. 


for 2h to stabilize the film resistivity. The cavity-patterned sili- 
con wafer was then anodically bonded with the patterned glass 
wafer, entrapping air at the required pressure. The bonded insula- 
tion wafers were then diced into individual chips of the required 
sizes. Fig. 8 shows the fabricated components of an integrated 
SR. 


3.4. Incorporation of commercial catalyst and assembly of 
an integrated SR 


The commercial BASF catalyst is incorporated in the form of 
a packed bed in the silicon microreactor bonded with glass. The 


(b) 


catalyst was received in a pellet form, and was ball milled and 
sieved to obtain catalyst particles of mean diameter 70-80 um. 
These particles were then introduced into the reaction zone 
through the inlet via while creating suction at the outlet. The 
reactor has an array of catalyst filter structures at the outlet with a 
gap spacing of 30 um ensuring that catalyst particles larger than 
this remain trapped within the reaction zone. Fig. 8(c) shows the 
silicon microreactor packed with catalyst. 

Next, the insulation chips are bonded on the front and back 
sides of the catalyst-loaded microreactor using high temperature 
epoxy. Fig. 9 shows the microreactor bonded with the insulation 
chips. 


Fig. 9. SR microreactor bonded with the insulation chips. (a) Frontside of SR microreactor bonded with insulation chips. (b) Backside of reactor bonded with 


insulation chips. (c) Assembled integrated microreformer (six-chip stack). 
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4. Results and discussion 
4.1. Thermal characterization 


Thermal characterization of an integrated reformer was car- 
ried out without flows and reaction in order to measure critical 
thermal parameters including heat requirement, heat transfer 
rates between components, lateral and vertical temperature dis- 
tributions, insulation effectiveness, and heat loss mechanisms. 

Electrical connections were made by attaching thin wires to 
the contact pads with conductive silver epoxy. All thin-film tem- 
perature sensors were first calibrated from room temperature to 
300 °C to relate the change in resistance with temperature. A lin- 
ear trend was observed by plotting normalized resistance (R/Ro) 
versus the temperature. The temperature coefficient of resistance 
(TCR) was then calculated using a linear fit, which allowed 
estimating temperatures from the experimental resistance data. 

For thermal characterization, the integrated device was not 
supported by any fixture to avoid external thermal contact. 
Instead, the device was suspended in air only by the attached 
wires as shown in Fig. 10. 

For characterization, electrical heating with the integrated 
thin-film heater was carried out to introduce a well-defined rate 
of heating. Heating was carried out until the required 260 °C 
temperature was obtained in the SR reaction zone. Once the 
device reached steady state, the current and voltage across the 
heater were measured along with resistances of temperature sen- 
sors using digital multimeters. The power required to maintain 
the temperatures of the reaction zone and each component was 
thus obtained. We were able to safely introduce close to 6.8 W 
using this heater. At this power, however, it was only possi- 
ble to obtain a temperature of 230°C in the reforming reaction 
zone because of the high losses. Furthermore, we observed that 
increasing the power above ~6.8 W resulted in a short circuit 
of the heater resistor and cracks in the glass from the combined 
effect of high current density and high temperature of ~550 °C 
achieved within the small heater zone. Fig. 11 shows the tem- 


Fig. 10. Thermal characterization of an integrated device suspended in air. 
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Fig. 11. Temperature obtained in different components vs. the power input to 
the uninsulated integrated device suspended in the ambient. 


peratures of different components of the integrated device as a 
function of the power input to the heater. 

Despite sub-optimal SR temperature reached under these 
condition, it was possible to obtain the required thermal isolation 
between the combustor, the reformer and the PrOx units. For the 
reformer temperature of 230 °C, the temperature obtained in the 
combustor and the PrOx were 564 and 150°C, respectively. The 
integrated device was designed to attain 350 and 155 °C, in the 
combustor and PrOx respectively for the reformer temperature 
of 260°C. The reaction zone remained at uniform temperature 
as measured by three temperature sensors ((Tr2, Tr3, Tr4) placed 
at different location on the backside. This uniformity of temper- 
ature results from the high thermal conductivity of the silicon 
reactor embedded between insulation chips. 

However, for the free surfaces open to the ambient, large 
lateral thermal gradients were observed even from these 
small-scale devices. For example, we observed a temperature 
difference of greater than 100 °C between a combustor (Te, mea- 
sured by a thin-film sensor placed in the center of the combustor) 
and 7,1, thin-film temperature sensor that was placed at a lat- 
eral distance of 600 um from the combustor. Therefore, it was 
necessary to measure the surface temperature distributions at 
even higher resolution than originally planned. In addition to 
thin-film temperature sensors, manually placed thermocouples 
were used to obtain temperatures at locations intermediate to the 
micro RTD sensors. 

Various heat loss pathways were identified and the individual 
contribution of these pathways was determined from the mea- 
sured heat input as a function of surface temperature. Heat is 
lost from the free surfaces by free convection and radiation. In 
addition to this, heat loss through attached suspending wires and 
thermocouples was considered. Therefore, at steady state, 


Qin = VI = Oross = Owires + Oconv. F Qrad. (4) 


Heat loss from these wires on each surface was estimated by 
applying the so-called pin fin theory [41]. Heat dissipated from 
the wires was calculated by integrating the loss over the side 
surface of the wires, yielding 


pa (hP'(T, — Tao )tanhal)n Š 


a 
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where a = hP'/kAw, P’=2z7R, is the perimeter of the wire, 
R is the radius of the wire (250 um), Aw the cross-sectional 
area of the wire, A the heat transfer coefficient = 50 W m~? K7! 
(assumed), T; the edge temperature of the respective surface, 
Too the ambient temperature, /= 10 cm, the length of wire, and 
n is the number of wires on each free surface. 

Since the substrate media incorporated here are transparent 
to significant portions of the emission spectrum involved, the 
spectral volumetric emission and absorption in silicon and glass 
were taken into account. For an isothermal, absorbing and emit- 
ting medium, the intensity 7} (L) of a monochromatic beam of 
radiation is given by the following equation [42]. 


BCL) = hoje" + In, — eo") (6) 


The first term, J;,(0) eka relates to the intensity of an 
impinging beam and its reduction within the solid medium of 
thickness L due to absorption. The second term, Jp, (1 — e kab), 
is the volume-generated emission taking place over the thick- 
ness of material, L. For a composite of adjacent silicon and 
glass pieces of thickness L and L’, respectively, equation (6) 
then can be used to determine intensity at the free surfaces of 
the composite. Equation (7), 


Si Py 
DL’) = Ig — eR E24 K 15) (7) 


gives the intensity at the free surface of the Pyrex™ piece in 
a composite of silicon and Pyrex™. This equation can also be 
rewritten to obtain the intensity J, (L) from the free surface of the 
silicon piece. In these equations, J}, = Ep,/7, where Ep, is the 
i P 
spectral blackbody emissive power (W m~? m7!) [37]. kaid 
are the spectral absorption coefficients of silicon and Pyrex™, 
respectively. These calculations are made for the wavelength 
range from 0.1 to 100 um to obtain the total emissive power E 
(W m7’), which is then multiplied by area to obtain Qraq. After 
determination of the radiation rate, the estimated values of Qwires 
and Q;aq were then subtracted from the measured heat loss to 
evaluate Qcony, aS in equation (8). 


Oconv = Qloss m Q wires = Qrad (8) 


Fig. 12 shows the steady-state heat loss due to conduction 
through wires, natural convection, and radiation (surface radia- 
tion and volumetric emission) as a function of SR reaction zone 
temperature. The heat loss estimate shows that the surface loss 
by what we initially considered to be typical natural convec- 
tion makes up the major fraction of overall loss at atmospheric 
pressure. 

An unexpected finding from these results is that the effective 
convective losses were found to be significantly higher at this 
length scale than the predictions obtained by correlations for 
objects of conventional scale. We found that the heat transfer 
coefficient for natural convection increases dramatically as the 
surface dimensions decrease and was severely under predicted 
in our initial estimates. Based on additional experiments and 
analysis, we developed the following dimensionless correlation 
to predict the effective heat transfer coefficient of natural con- 
vection for micro to mesoscale objects (sub-mm to mm scale). 
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Fig. 12. Total steady-state heat loss measured experimentally and estimated heat 
loss by various pathways vs. the SR reaction zone temperatures. 


The detailed data analysis and derivation of this correlation can 
be found elsewhere. 


Nyu = 1.7Npa?®8 (9) 


where Nyu is Nusselt number and Nga Rayleigh number. The 
uncertainty in predicting Nusselt numbers from these correla- 
tions is approximately +15% attributed mainly to experimental 
error and assumptions made during estimation of the different 
losses. 

The observed high convective losses to the ambient are 
believed to be the combined effect of air conduction and con- 
vective transport with the former dominating at this scale. When 
we compared the pure convective heat rate using a conventional, 
macroscale correlation from the literature [37], we found that it 
could only account for 30% of the effective transport rate. This is 
in agreement with the results of others who found extraordinarily 
high effective rates of convection from mm- to cm-scale surfaces 
[43,44], leading to suggestions that direct conduction in the air 
is responsible for the discrepancy. We adhere to this hypothesis 
as well, and hence the correlation we developed (equation (9)) 
conveys an effective convection coefficient, since we believe 
another process besides convection must be contributing. Gas- 
phase (air) conduction could be explained as significant with 
decreasing system size in the following way: the small volume 
of gas heated by a small (mm- to cm-scale) interfacial area is 
changed in density due to an increase in temperature. While 
with macroscale objects, this density difference drives the flows 
associated with natural convection, the buoyancy forces on the 
small volumes of heated gas in this case are too small to strongly 
overcome the forces of viscous drag which oppose the convec- 
tive motion. As a result, the gas does not flow significantly, and 
direct conduction from the surface becomes dominant. 

Thermal radiation is also important in the consideration of 
our measurements as it amounts to a significant fraction of loss 
especially at high temperatures. However, from experimental 
measurements of steady-state heat loss and derived correlations, 
we observe that gas conduction and convection dominate loss by 
radiation even at fairly high temperatures for small-scale devices 
compared to macroscale devices. For all temperature ranges 
measured, heat loss via the attached wires was less than 7%. 
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Fig. 13. Heat loss by air conduction and natural convection as a function of 
the SR reaction zone temperature for the integrated device suspended in air and 
under vacuum. 


Thermal characterization of this assembled device was then 
carried out under vacuum to determine the effect of sub- 
atmospheric pressure on steady state heat loss. The integrated 
device resting on the attached wires was placed in the vaccum 
chamber. The chamber was evacuated down to 50 mTorr, and 
the steady-state heat losses were characterized as a function of 
surface temperatures. 

Measurements under vacuum showed that ~4.3W was 
needed to obtain the required 260°C temperatures in the SR 
reaction zone. The steady state heat loss data was then analyzed 
to determine relative magnitudes of conduction through wires, 
air conduction plus natural convection, and radiation losses. 
Fig. 13 shows the comparison of heat losses by air conduc- 
tion plus natural convection as a function of SR reaction zone 
temperature for the ambient pressure and vacuum tests. For the 
SR reaction zone temperature of 230°C, the heat loss by air 
conduction plus natural convection was 2.1 W under vacuum 
compared to 4.8 W observed at ambient pressure. Therefore, it 
was possible to reduce this loss by about 60%. 

Additional thermal characterization experiments were car- 
ried out both at atmospheric pressure and under vacuum where 
the integrated device was covered in commercial fiberglass insu- 
lation. Table 3 shows the heat input and temperatures attained 
in combustor, SR, and PrOx zones, and a breakdown of steady- 
state heat losses into convective and radiative losses for each 
experiment. 


Table 3 


189 


As shown in Table 3, the device when covered in loosely 
packed fiberglass insulation showed the best performance at 
both atmospheric and vacuum tests. Losses decrease with an 
increase in thickness of the insulation, resulting from reduced 
temperature gradients and heat transfer coefficient. However, 
using a large quantity of insulation adds to the weight and espe- 
cially volume of the overall system, which translate into reduced 
energy density. 


4.2. Reaction characterization 


Experiments are carried out to determine conversion and 
selectivity, hydrogen production rate (corresponding power out- 
put), heat required to drive the reforming reaction, and thermal 
effects on reaction. 

The integrated microreformer was characterized on a custom- 
built test setup. A high performance liquid chromatography 
(HPLC) pump was used to introduce the premixed liquid reac- 
tants of water and methanol to the vaporizer. The vaporizer was 
built by wrapping flexible heaters around the stainless steel tube 
of 1/16 in outside diameter. The vaporized mixture was then fed 
to an inlet of the silicon microreformer placed on the interface 
block. Separate inlet lines are provided to introduce measured 
flows of nitrogen and 5% hydrogen to the microreformer for the 
catalyst reduction. The stainless steel reactor block was designed 
to ensure minimum thermal contact with the device. The prod- 
uct was sent to the condenser to remove unconverted water and 
methanol. A downstream mass flow meter was used to measure 
the flow rate of the dried reformate stream before it was analyzed 
by micro gas chromatography (Varian CP-4900). The whole 
tubing system was wrapped with tape heaters and insulated to 
prevent condensation. 

Electrical connections to the exposed contact pads on the inte- 
grated thin-film heater and temperature sensors were made from 
outside using standard alligator clips. The vaporizer heaters, 
mass flow controllers, and the GC were automatically monitored 
and controlled by synchronizing the setup with a computer using 
LabView software. 

Prior to the reaction, the catalyst was reduced with 2 sccm of 
5% H2/N2 at 280°C for 2h. The reforming reaction was then 
carried out by introducing a methanol-water mixture with a 1:1.3 
molar ratio at a defined rate along with a stream of 5 sccm of 
N2 as an inert chromatography reference. The flow rates of H2, 


Power required to maintain reformer temperature at 260°C, corresponding to temperatures obtained in the combustor, SR, and PrOx zones, and breakdown of 


steady-state heat loss into convective and radiative losses for each experiment 


Experiment Qin (W) Tcomb CO) Tsp CC) Tprox CC) Qeon. (W) Qrad (W) 
Suspended on wires in ambient 6.8 564 230 149 4.8 1.7 
Suspended on wires under vacuum 4.3 452 260 176 2.4 1.7 
Covered with fiberglass insulation of size, /=4 cm, w=2.3cm, t=2.4cm, 5.1 532 260 166 3.6 1.3 
in ambient 
Covered with fiberglass insulation of size, /=4 cm, w = 2.3 cm, t=2.4cm, 2.8 414 260 186 1.3 1.3 
under vacuum 
Placed in insulation of size, l= 10 cm, w = 8 cm, t= 8 cm, in ambient 2.5 411 260 181 1.3 1.0 
Placed in 10cm x 8cm x 8 cm insulation, under vacuum 1.8 386 260 181 0.6 1.0 


The convective loss shown here actually accounts for the combined effect of air conduction and natural convection losses. 


190 K. Shah, R.S. Besser / Journal of Power Sources 166 (2007) 177-193 


100,00 
90,00 
80.00 
70.00 
60.00 
50.00 
40.00 
30.00 
20.00 
10.00 

0.00 


Methanol Conversion, % 


0 20 40 60 80 «100 120 140 160 180 200 
Time, min. 


Fig. 14. Methanol conversion obtained at different temperatures. Methanol- 
water feed: 0.009 ml min™!, S:C = 1.3:1, WHSV: 5.1 h7}, GHSV: 11670h7!. 


CO7, and CO were obtained by relating their GC intensities to 
the signal from the 5 sccm N2 flow. 

Because of the relatively large thermal mass of the block and 
high losses, it was not possible to heat the device to the required 
temperature with the device’s integrated thin film heater alone. 
Hence a tape heater was used to maintain the test block at a tem- 
perature of 150-200 °C. The integrated microreformer was then 
heated to a desired temperature in the SR reaction zone using the 
thin-film heater. The reaction was carried out at different tem- 
peratures from 220 to 280°C with a combined feed flow rate 
of 0.009 ml min~! (methanol flow rate of 1.36e—4 mol min~!). 
Fig. 14 shows the methanol conversion obtained at different 
temperatures. At 280°C, close to 100% methanol conversion 
was achieved, producing 9 sccm of H2. At the reaction temper- 
ature of 220°C, the methanol conversion decreased to 50%, 
reducing the H2 produced to 5.2sccm. CO was not formed 
in detectable amounts at the reaction temperature of 220°C, 
however, it increased to 1.1% as the temperature increased to 
280°C. This CO can be further minimized to less than 10 ppm 
by employing PrOx microreactor to prevent poisoning of fuel 
cell electrodes. 

We did not observe any major change in temperature pro- 
files during the reaction arising from either the flows or the 
endothermic nature of the reaction. This occurred because the 
heat required to carry out the endothermic SR reaction is quite 
small relative to the amount of heat to be input (because of high 
heat losses) in order to maintain the required reaction tempera- 
ture. 

The effect of feed flow rate (and equivalently space veloc- 
ity) on the methanol conversion was determined by changing 
combined feed flows from 0.007 to 0.015 ml min7! at a reac- 
tion temperature of 260 °C. Methanol conversion was decreased, 
while the reformate and hydrogen flow rates increased with 
increasing feed flow rate (resulting in decreasing residence 
time). With the total feed rate of 0.015mlmin~!, 82% of 
methanol conversion was achieved, producing 14.3 sccm of 
hydrogen. Considering 60% FC efficiency and 75% H; utiliza- 
tion, this stream is sufficient to produce 1.04 W from the FC. 
Fig. 15 shows the flow rates of H2, CO2, and CO in the refor- 
mate stream as a function of feed flow rates. Also shown are the 
corresponding anticipated powers from the FC. After the sepa- 
ration of condensibles, the concentrations of H3 in the reformate 


Hy: 143 sccm 
1.04 W = 


Hy: 11.6 sccm 


Ha: 6.7 sccm 


o5 wW 


Flow rates of the reformate species, 8cem 


0.007 0.009 0.012 0.015 
Methanol-water feed flows, ml min-1 


Fig. 15. Flow rates of H2, CO2, and CO produced as a function of methanol- 
water feed rate. 


stream remained constant at 74% at all the conditions of tem- 
perature and feed flow rate. The catalyst remained stable over 
a reaction period of 6h studied here, however, reaction should 
be carried out over greater lengths of time to evaluate catalyst 
performance and degradation issues. 

Microreactors are known to have faster heat and mass trans- 
fer rates than conventional fixed bed reactors and are normally 
assumed to operate in the reaction-limited mode [45]. How- 
ever, for the practical reactor system considered, it is necessary 
to check the influence of mass transfer on the reaction perfor- 
mance. We assumed the absence of internal mass transfer as 
it can be safely neglected for catalyst particles with a diame- 
ter smaller than 100 um [36]. External mass transfer, involving 
mass transport of reactants from the bulk to the catalyst surface 
usually becomes important in microreactors. Based on the calcu- 
lations of the diffusion coefficients, it was shown that methanol 
is most likely the limiting reactant in the SR reaction [46]. 
We employed Mears’ criterion to examine the methanol mass 
transport limitation. 

According to Mears, there is no mass transfer limitation if 
the following criterion is fulfilled [47]: 


d 
MORS. 3 (5 (10) 
2k gCMeoH 


where rmeou is the measured reaction rate, dp the hydraulic 
diameter of the reaction zone, kg the mass transfer coefficient, 
and rmeou is the methanol concentration in the feed. 

Due to the small hydraulic diameter and low flow rates 
involved, the flow in the microreactor was laminar. The mass 
transfer coefficient was determined by considering that in 
the fully developed laminar flow, the Sherwood number Sh 
approaches a constant value. 

Sh= Kedr = 3.66 (11) 
DAB 

where Dap is the diffusivity of methanol in the gas mixture, 

calculated on the basis of mole fractions of the species in the 

gas mixture [48]. 
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Fig. 16. Methanol conversion predicted by simulation and measured experimentally (a) as a function of feed flow rate with a reaction temperature of 260 °C, (b) as 


a function of temperatures with a feed rate of 0.009 ml min~!. 


The value of rmeou dn/2kgCMeoH estimated was 7.5e—2, 
which implies that there was no external mass transport lim- 
itation. Fig. 16(a) and (b) shows the comparison between the 
methanol conversion measured and predicted by the kinetic 
model at different temperatures and flow rates, respectively. 
The experimental conversion obtained matched reasonably with 
that predicted by a kinetic model based on the literature kinetic 
expression, neglecting mass transfer resistance. This further 
corroborates the absence of external mass transport limitation, 
confirming that the reaction was rate limiting. 


4.3. Packaging suggestions and scale up 


Based on the understanding gained from thermal characteri- 
zation experiments, packaging suggestions are made to reduce 
air conduction, convection, and radiation losses. 


- In scaling up, it is desirable to form a stacked three- 
dimensional assembly in preference to using larger layers (3D 
architecture instead of planar sheets). This reduces the surface- 
to-volume ratio and thus the surface convective and radiative 
losses. 

- Heat losses by air conduction and natural convection consti- 
tute a significant fraction of total heat loss at ambient pressure. 
These heat losses were reduced by covering the device in insu- 
lation and operating under vacuum. Use of low-k insulating 
materials helps reduce these losses. However, there is a trade- 
off on the amount of insulation that can be used without adding 
much to the weight and volume of the integrated system. Addi- 
tionally, it would be helpful to use a low-k gas like argon 
instead of air surrounding the integrated device to reduce gas 
conduction losses. We observed that the use of vacuum below 
50 mTorr reduces these losses by at least 60%. However, at this 
level of vacuum, the air conduction and convection could not 
be entirely suppressed. This is consistent with other reports 
in literature [32,35]. Vacuum packaged insulation is superior 
for reducing these heat losses without significantly increasing 
volume. However, it seems that a vacuum down to 10~* Torr 
or below is required to minimize losses through the ambient 
air [32]. This presents a challenge in developing a reliable 


and hermetic vacuum package that can hold such a low pres- 
sure for a long period of time. We suggest employing both 
low-k insulating material and vacuum packaging to reduce air 
conduction and convective losses. 

- Radiative heat losses from a free surface can be reduced by 
using reflective insulation or by employing high reflectivity 
and low emissivity materials. 


We observed that heat loss to the ambient via gas conduc- 
tion, convection, and radiation was substantial compared to 
the required heat for critical energy-consuming FP functions 
including fuel vaporization and the endothermic SR reaction. 
Considering the integrated device of the vaporizer, combus- 
tor, SR, and PrOx having a size of 1cm x 1cm x 0.5cm, if 
we employ the above packaging recommendations, it would 
still result in a minimum heat loss of 1.0 W. The total energy 
required considering heat loss (1.0 W) plus energy to drive the 
endothermic SR reaction (0.13 W) and fuel and water vapor- 
ization (0.20 W) is then 1.33 W. Considering 75% hydrogen 
utilization in the FC, hydrogen recycled from the FC anode when 
burned gives only 0.38 W, thus an additional 1.33-0.38 = 0.95 W 
would required by combusting methanol fuel. The device can 
produce 3.84e—4 mol min! of hydrogen, resulting in a lower 


Total heat required (heat loss + heat of reaction + heat 
required for sensible heating and vaporization of 
MeOH-water feed), W 


LHV of un-utilized 
hydrogen, W 


Heat, W 


Anticipated heat loss. W 


0 - 
0 5 10 15 20 
Fuel cell power output, W 


Fig. 17. Anticipated heat loss, reformer energy requirement, and LHV of the 
un-utilized hydrogen in the anode-off gas as a function of FC power output. 
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heating value (LHV) of 1.55 W, while considering FC efficiency 
as mentioned previously it can produce a net 0.65 W output. Thus 
we will begin to see the benefits of an integrated FP/FC approach 
as systems are designed for powers exceeding a few W. Fig. 17 
shows the anticipated heat loss after stacking of reactor layers 
and employing packaging recommendations, including the heat 
duty required for endothermic SR reaction and for fuel and water 
vaporization, and the LHV of un-utilized hydrogen available at 
different power outputs of 1, 10, and 20 W from the FC. 


5. Conclusions 


Understanding thermal management issues in a portable FP 
is crucial in order to reduce heat losses and improve thermal 
efficiency. We implemented silicon microchemical systems to 
investigate thermal integration of FP reactor components with a 
methanol SR as a model. The fabricated model components of 
an integrated microchemical device were used for the acquisi- 
tion of critical thermal parameters of the overall coefficients for 
the transfer of heat from components to one another and to the 
outside environment by various heat loss pathways. Based on 
thermal characterization experiments, the heat loss mechanisms 
and the effective convective heat coefficient were determined 
and suggestions were made for scale up and implementation 
of packaging schemes to reduce different modes of heat loss. 
Acquisition of the thermal experimental data helped to ver- 
ify and improve the model used to simulate the system. This 
combined experimental/modeling approach can then be used 
as the quantitative basis for designing an integrated portable 
FP. 

The SR microreactor produced 14.3 sccm of hydrogen (cor- 
responding to a net power output of 1.04 W from the FC) with 
relatively low levels of carbon monoxide. The results effectively 
demonstrate the potential of microreactor-based on-demand H2 
generation, and provide impetus for a subsequent phase of devel- 
opment involving completely integrated subsystems. 
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